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We propose to use a laser to generate magnetic Bloch oscillations in one dimensional easy-axis 
ferromagnets at low temperatures. This proposal is investigated numerically in details for material 
parameters relevant for C0CI2 ■ 2H2O. 



PACS numbers: 75.78.-n,75.40.Mg,75.10.Pq,78.47.D- 

According to quantum mechanics a particle in a pe- 
riodic potential will oscillate in response to a constant 
force. Such Bloch oscillations (BO) were predicted 
in the early days of quantum mechanicsii^, but have 
only recently been experimentally demonstrated in very 
clean semiconductor superlattices^ and in Bose-Einstein 
condensates'*. 

In condensed matter systems the particles need not be 
of the ordinary kind resembling electrons. In particular, 
an elementary particle in a one dimensional anisotropic 
ferromagnet is a domain-wall separating regions of up 
and down spins. Such domain- walls can have a dispersion 
relation like that derived from a periodic potential, and in 
the presence of a uniform magnetic field these magnetic 
systems have been predicted to show BO^. In particular 
the blue crystalline material C0CI2 • 2II2O and CoNb2 06 
have been proposed as candidate materials. However, no 
BO have been experimentally observed in these to date. 

In a recent article*' we have revisited the material 
C0CI2 • 2H2O and refined the predictions of Ref. [l for 
detecting spectral signatures of BO in neutron scattering 
experiments, taking also into account extra interactions 
present in the material. While our result indicates that 
the spectral signatures of BO can indeed be observed in 
neutron scattering, the signatures are relatively weak, at 
the 10% level of the total spectral weight at finite tem- 
peratures. This is a consequence of neutrons being a 
relatively weak probe as they cause only single spin-flip 
excitations. 

We propose here a more direct way to generate BO 
by keeping the material at low temperature and in- 
duce excitations using a short laser-pulse. Upon turn- 
ing off the laser-pulse the magnetization of the mate- 
rial will continue oscillating at the Bloch frequency. It 
has been known since long ago2*S that light in the far- 
infrared frequency range can induce magnetic excitations 
in C0CI2 • 2H2O, but no time-dependence of the magne- 
tization was studied there. 

In this article we model the laser-pulse as a time- 
dependent perturbation to the Hamiltonian and investi- 
gate its effects by solving the time-dependent Schrodinger 
equation numerically. We show that BO can be gener- 
ated this way, and give appropriate laser frequencies and 
pulse-duration times. 

The magnetic properties of C0CI2 -2112 are described 



by the spin- 1/2 Hamiltonian 

^ = - E (^-^'^^+1 + h-St) + Ha, (1) 

i 

where denotes subdominant terms to be discussed 
below. The ferromagnetic coupling — 36.5 is the 
dominant term in the Hamiltonian. Alone it causes 
neighboring spins to align their spin z-components, thus 
the energy of an excited state depends on the number 
of anti-aligned spin neighbors; domain-walls, where each 
domain- wall costs an energy oi Jz/'i-- In the presence 
of an external magnetic field along the z-axis, h^, the 
energy will also depend on the number of spins oppos- 
ing the field, implying pair-wise confinement of domain- 
walls. Such a bound-state of two domain-walls separating 
I overturned spins, known as a spin cluster excitation^, 
or simply a domain, has an energy + h^l above the 
ground state energy. 

The term Hd describes additional couplings that partly 
give dynamics to the domain walls and partly induce 
more domain- walls, 

i 

+J^{StSr^, + SrSl,)], (2) 

where Ja = 3.8 and J_l = 5.43 K^. 5*^ = Sf ±iSl are 
the usual raising and lowering operators. The Ja-term 
can move a domain- wall two lattice spacings, thus mixing 
states with even (or odd) I. Similarly Jj_ gives kinetic 
energy to the I = 1 domain state. Both of these terms 
can also induce new domain- wall pairs. However, with 
being the dominant coupling, extra domain- walls will be 
energetically costly, thus we will restrict our calculations 
to states having a small number of domain- walls Ndw ■ 

When restricting to Ndw < 2, iJ can be diagonalized. 
The energy spectrum is En — Jz + IJ^nhz, where /l(„ is 
found by solving an equation involving a ratio of Bessel 
functions^. For high energies /i„ ~ n, an integer, thus the 
energy spectrum becomes equidistant. The correspond- 
ing energy eigenfunctions {ipn) are Bessel functions^. 
From these one can construct a time-dependent state 
lx(^)) — J2n Qne~*'^"*|V'n)- When this state is dominated 
by energy eigenstates with energies from the equidis- 
tant region, the time-dependence of the magnetization 
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FIG. 1: In C0CI2 • 2H2O the spin z-axis corresponds to the 
crystallographic b axis. The strong ferromagnetic Jz couples 
spins along the c axis. A laser beam is shown propagating at 
normal incidence to the b-c-plane with magnetic field polar- 
ization in the b direction. 



MM = {x{t)\Y.;Snx{t)) becomes 

Mz{t) = c - 2a;o ^ i?e (a*a„+2e 
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where c is a constant, ujb = ^h^ is the Bloch frequency, 
and xq = 2Ja/hz- Thus the BO amplitude is propor- 
tional to 4:Ja/hz times a factor which depends on the 
probability amplitudes a„ of the excited states. Our BO 
amplitude is a factor 2 larger than obtained in Ref. ^ as 
we consider the size oscillations of a single domain having 
two domain- walls, while they considered the motion of a 
single domain- wall. 

In order to populate the excited levels at very low tem- 
peratures, thus producing BO, we propose here to use a 
laser with a wavelength in the far-infrared. The Co elec- 
trons causing the magnetism are d-shell electrons thus 
having no electric dipole moment. We therefore model 
the laser as an extra time-dependent magnetic field which 
couples to the spins as 



(4) 



where ut is the laser frequency and Bq is the laser mag- 
netic field amplitude. We have here assumed a linearly 
polarized laser beam such that the magnetic field is along 
the Ising direction. This corresponds to the crystallo- 
graphic b direction in C0CI2 • 2H2O. Such a setup can be 
made by cleaving the crystals in the b-c plane and direct- 
ing the laser at normal incidence to this surface polarizing 
the laser beam such that the magnetic field points along 
the b-direction, see Fig. [T] We have assumed the laser 
beam to be coherent along its front, and also through the 
crystal. The laser considered here have a wavelength of 
about 0.3 mm thus for this approximation to be good the 
crystal should be thinner than this. For thicker crystals 
there will be an additional phase-shift associated with 
the depth. 

The time-dependent perturbation can be treated nu- 
merically for large system sizes when restricting to states 
where Ndw < 2. This restriction implies that the energy 



gap between the ground state and any excited state will 
depend on the system size iV. This can be understood 
by considering the perturbative energy correction from 
virtual processes involving the creation and destruction 
of an additional domain. As there are roughly N places 
to insert the new domain, the energy correction will be 
proportional to N. When restricting to Ndw < 2 the 
ferromagnetic state receives this correction, but not the 
states having one domain, as their corrections come from 
the excluded Ndw > 4 sector. In order to make the en- 
ergy gap intensive we redefine the coupling between the 
ferromagnetic and the 1 — 2, Ndw = 2 state in the Hamil- 
tonian by dividing it by a factor ViV jS, This effectively 
makes the correction to the ferromagnetic state indepen- 
dent of system size. 

Starting in the ground state of H, the time-dependent 
Schrodinger equation is solved iteratively numerically 
with the laser field He^t present. In the iterations we keep 
the 300 lowest energy states of H with zero momentum 
and Ndw < 2. The laser frequency to is tuned such that 
Lo = En—Eo where n corresponds to an energy level in the 
region where the spectrum is approximately equidistant. 
We choose n — 12 corresponding to E12 — i?o ~ l-BJz f^or 
a static magnetic field hz — 0.05 J^. In practice when us- 
ing a laser with a fixed wavelength, resonance can instead 
be found by changing the static magnetic field thereby 
adjusting the energy levels. The iterative solution gives 
time-dependent amplitudes of the different energy levels. 
Fig. [2] shows the probabilities of finding the system in se- 
lected levels as a function of time. Only even n states are 
excited because the S'^-terms do not flip any spins. The 
black dot-dashed curve in Fig. [2] shows how the ground 
state is depleted. The minimum of the ground state pop- 
ulation coincides with the maximum of the population 
of level n = 12, green solid curve, and occurs at a time 
T = TTh/ujR where ojr = A/|Bgai2p + (uj - {E12 - Eq)Y, 
where ai2 is the matrix element of Sf between the 
ground state and the n = 12 excited state. These oscilla- 
tions are in essence Rabi oscillations. Exciting the level 
71 = 12 alone does not give appreciable amplitude for 
Bloch oscillations as one also needs to populate the lev- 
els with n+2 (or n— 2), see Eq.|ni This can be achieved by 
using a large amplitude of the laser, we used Bq = 0.2 Jz, 
thereby causing off-resonant tunneling between the n and 
the n ±2 levels, see Fig. [2l red dotted and blue dashed 
curves. These off-resonant tunneling processes are fast, 
thus the population of the nearby levels follows closely 
that of the central level. 

In order to seek the maximum amplitude of BO we 
turn off Hext at the maximum population of the central 
excited level. A close look at the time dependence of 
a* a„+2 near the cutoff time, reveals that the dominating 
terms have the same phase, thus the amplitude of BO is 
proportional to ^„ |a*a„+2| which is shown in the lower 
inset of Fig. [2l From this we see that it oscillates fast 
with a frequency corresponding to E12 — Eq. It may 
be difficult to turn off the laser when this quantitiy is 
maximal. However, this is not a major concern as the 
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FIG. 2: (color online) Population |ajp of selected energy lev- 
els j, indicated in the legend, as a function of time after turn- 
ing on the laser, H^xt- = (^'12 — Eo)/h and Bq — 0.2 ■ 
The upper right panel shows a zoom in on the boxed time 
region. The lower right panel shows the time dependence of 
the sum of products of two nearby probability amplitudes, 
X^n l^^nOiii+al- The time averaged value is shown as the hor- 
izontal dashed line, t is measured in units of h/Jz which is 
0.2- 10"^^s for C0CI2 • 2H2O. 



time averaged value is about 75% of the maximum value. 

Turning off the laser at a time t when |a*a„-|-2| is 
maximal, and letting the system evolve further in time 
without Hext, produces the BO shown in Fig. [3l As our 
simulation only allows single domain excitations we have 
plotted the relative size of the domain, measured by the 
expectation value of the number of spins opposing the 
fi_eld Nii{t) = M^{t) - N/2 divided by its time average 
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n for excitation £"„. We see that the relative size 



of the domain oscillates between 0.6 and 1.4 correspond- 
ing to a size between 7.2 and 16.8 for N12 = 12. Thus 
the amplitude is 4.8 which is close to the expected value 
'iJa/hz X 0.6 — 5 from Eq. [31 Allowing a finite density p 
of coherently oscillating domain states, the relative size 
of a single domain state shown in Fig. [3] will be propor- 
tional to the experimentally relevant quantity, the time- 
dependent relative magnetization: {Mz{t) — M)/M = 
(p7Vu/(l/2 - pN,^){N,^{t)/N,^ - 1). 

The BO in Fig. [3] are not simple harmonic. A Fourier 
transformation of the beating pattern is shown in the 
lower panel of Fig. [S] Two peaks are clearly visible. 
They correspond to cji = E14 — E12 = O.lJz and 
UJ2 = E12 — Eio = O.llJ^. Thus the beating pattern 
is due to the deviation from an equidistant ladder spec- 
trum. This frequency difference can be made smaller by 
exciting higher energy bound states where the spectrum 
is closer to being equidistant. 

In this setup we rely on off-resonant tunneling in order 
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FIG. 3: Size oscillations of a domain excitation, measured as 
the number of down spins Nii/Nii vs. time after the laser 
is switched off (upper panel) , and its Fourier spectrum (lower 
panel) . 



to populate nearby levels. This requires a large laser 
amplitude. Instead one might use two small-amplitude 
lasers each in resonance with nearby levels. One can also 
change the laser polarization to have a component along 
the spin-x direction. This will induce transitions between 
the even and odd n states. Our simulations show BO also 
in this case, but now with more frequency components 
due to the deviation from equidistant spectra, as in this 
case both the even and the odd sectors participate. 

The restriction to Nij^ < 2 and the associated redefini- 
tion of the coupling to the ferromangetic state can raise 
doubts about the validity of the matrix elements calcu- 
lation, also it does not allow for any discussion of in- 
teractions between domains. We have therefore numeri- 
cally also investigated cases where we allow more domains 
without any redefinition of couplings. Computer perfor- 
mance restrictions let us consider Ndw < 6 for TV < 34. In 
Fig. m it is seen that the N dependence of the energy gap 
A to the first excited state decreases as higher domain 
wall sectors are included. For comparison the Ndw < 2 
result with redefined coupling to the ferromagnetic state 
is shown as the dashed line. 

In order to identify the elementary domain excitations 
and the associated transition matrix element of the laser 
field between these and the ground state, we construct 
the approximate single domain creation operator— 



-'^^'/V(z-,.„)/2(^ocosp)n£!,ri5+ (5) 



with momentum p — and let it act on the ferromagnetic 
state. The sum is restricted to even values of I and Jm is 
the Bessel function of the first kind of order m. We then 
compute the overlap of this with the exact eigenstates 
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FIG. 4: Energy gap A to the first excited state vs. system 
size for different Ndw The legend specifies the maximum 
Ndw included in the diagonalization. The dashed line shows 
the energy gap for Ndw < 2 with the l/\/iV redefinition of 
the coupling to the ferromagnetic state. 




FIG. 5: a) Excitation energy of elementary domain states vs. 
n (triangles). Also shown is \an\ vs. n (circles) on a semi- log 
scale. Filled symbols refer to Ndw < 6 and A'' = 34, while 
open symbols are for Ndw < 2 with redefined coupling to 
the ferromagnetic state, b) Interaction energy of two n = 12 
excitations vs. separation distance r for Ndw < 6 and A'' = 34. 



of the A^diD ^ 6 system. For each value of n we pick 
the state with maximum overlap. Fig. [S^) shows the 
energies of these states. We see that for high enough n the 
energies become equidistant and agree well with what we 
found for Ndw < 2. The laser transition matrix element 



between these states and the ground state behaves as 
VNan- The coefficients a„ are very close to those we 
found for Ndw < 2, see Fig.[5K), thus the Rabi frequency 
is not changed provided interactions between domains are 
negligible. Note that a„ drops very fast with increasing 
n. 

Insights about interactions can be gotten by identify- 
ing two-domain excitations. We construct approximate 
two-domain states as pairs with total momentum zero of 
two n = 12 single domain states separated by a distance 
r, with creation operator 6j „ = J2p (^p,n(^-p,n^^^^ ■ We let 
this act on the ferromagnetic state, retain only the terms 
having four domain- walls, and compute its overlap with 
the exact eigenstates of the Ndw < 6 system. The energy 
of the states with maximum overlap minus two times 
the single-particle excitation energy E12 are shown in 
Fig-Eb) as a function of r. We interpret this as the inter- 
action energy of domains separated by r lattice spacings. 
We find that the functional form a — b[l/r+l/{N — r)] fits 
the results reasonably with a = 0.227^ and b = 1.37 Jz- 
The positive a is caused by the restriction on Ndw which 
tends to overestimate energies in higher domain- wall sec- 
tors relative to those in lower sectors. This interaction 
energy causes an inhomogeneous broadening of the reso- 
nance frequency, and can lead to an upper limit on the 
density of domains excited by the laser together with 
an increased Rabi frequency due to interaction blocking 
cffectsiSiii. It can also change the BO frequency locally. 
However, in order to study this one needs to look at dif- 
ferences between interaction energies for neighboring en- 
ergy levels. We have not been able to study long enough 
chains to address this issue with sufficient precision. The 
lifetime of a domain excitation can also be reduced due 
to collisions, however the domain excitations are heavy 
due to their flat dispersion, so we expect a significant 
lifetime reduction only at a high density of excitations 
when neighboring domains are touching. 

We conclude that it should be possible to excite mag- 
netic BO in C0CI2 • 2H2O using a laser at low tempera- 
tures in a static magnetic field. In addition our simula- 
tions using the Hamiltonian and parameters from Ref. [12 
also indicate that BO in CoNb206 may be generated in 
a similar way. 
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